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Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and 
their associated instabilities (e.g., Buneman, Weibel and other two-stream instabilities) created in 
collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 
3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration 
associated with a relativistic jet front propagating into an ambient plasma. We End that the growth 
times of the Weibel instability in electron-positron jets are not affected by the (electron-positron 
or electron-ion) ambient plasmas. However, the amplitudes of generated local magnetic fields in 
the electron-ion ambient plasma are significantly larger than those in the electron-positron ambient 
plasma. The small scale magnetic field structure generated by the Weibel instability is appropriate 
to the generation of “jitter” radiation from deflected electrons (positrons) as opposed to synchrotron 
radiation. The jitter radiation resulting from small angle deflected electrons may be important for 
understanding the complex time structure and spectral evolution observed in gamma-ray bursts and 
other astrophysical sources containing relativistic jets and relativistic collisionless shocks. 


I. INTRODUCTION 

Nonthermal radiation observed from astrophysi¬ 
cal systems containing relativistic jets and shocks, 
e.g., active galactic nuclei (AGNs), gamma-ray bursts 
(GRBs), and Galactic microquasar systems usually 
has power-law emission spectra. In most of these 
systems, the emission is thought to be generated by 
accelerated electrons through the synchrotron and/or 
inverse Compton mechanisms. Radiation from these 
systems is observed in the radio through the gamma- 
ray region. Radiation in optical and higher frequencies 
typically requires particle re-acceleration in order to 
counter radiative losses. 

Fermi acceleration is the mechanism usually as¬ 
sumed for the acceleration of particles in astrophysical 


environments characterized by a power-law spectrum. 
This mechanism for particle acceleration relies on the 
shock jump conditions in relativistic shocks [e.g., 1, 
2]. Most astrophysical shocks are collisionless since 
dissipation is dominated by wave-particle interactions 
rather than particle-particle collisions. Diffusive shock 
acceleration (DSA) relies on repeated scattering of 
charged particles by magnetic irregularities (Alfven 
waves) to confine the particles near the shocks. How¬ 
ever, particle acceleration near relativistic shocks can¬ 
not be characterized as DSA because the propaga¬ 
tion of accelerated particles ahead of the shock can¬ 
not be described by spatial diffusion. Anisotropies in 
the angular distribution of the accelerated particles 
are large, and the diffusion approximation for spatial 
transport does not apply [3]. 

Particle-in-cell (PIC) simulations can shed light on 
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FIG. 1: 2D images in the x — z plane at y = 43A ((a) and (b)) and in the x — y plane at 2 = 230A a flat jet injected into a 
unmagnetized ambient medium shown at t = 28.8/u} pe . Colors indicate the y-component of the current density (J y ) ((a) 
and (b)) and z-component of the current density (J z ) ((c) and (d)) for electron-positron ((a) and (c)) and electron-ion 
((b) and (d)) ambient plasmas. Their peaks are (a) 15.6, (b) 23.9, (c) 54.7, and (d) 95.6. Arrows show J z , J x ((a) and 
(b)) and B x , B y ((c) and (d)). 


the physical mechanism of particle acceleration that 
occurs in the complicated dynamics within relativistic 
shocks. Recent PIC simulations using injected rela¬ 
tivistic electron-ion jets show that acceleration occurs 
within the downstream jet, rather than by the scat¬ 


tering of particles back and forth across the shock as 
in Fermi acceleration [4-10]. In general, these inde¬ 
pendent simulations have confirmed that relativistic 
jets excite the Weibel instability [11]. The Weibel in¬ 
stability generates current filaments with associated 
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FIG. 2: One-diinensional cuts along the z-direction (25 < z/A < 154) of a flat jet. Shown are the ai-components of 
the magnetic field shown at t = 238.8/a; pe for unmagnetized electron-positron (a) and electron-ion (b) ambient plasmas. 
Cuts are taken at x/A = 38 and y/A = 33 (blue), 43 (red), 53 {green) and separated by about an electron skin depth. 


magnetic fields [12], and accelerates electrons [4-10]. 

In this paper we present new simulation results of 
particle acceleration and magnetic field generation in 
relativistic electron-positron shocks using 3-D rela¬ 
tivistic electromagnetic particle-in-cell (REMP) simu¬ 
lations. We have performed two simulations in which 
an electro-positron jet is injected into two different 
ambient plasmas (electron-positron and electron-ion). 
In our new simulations, the growth rate of the Weibel 
instability and its evolution with different ambient 
plasmas have been studied without an initial ambi¬ 
ent magnetic field. 


II. SIMULATION SETUP AND RESULTS 

Two simulations were performed using an 85 x 
85 x 320 grid with a total of 180 million particles 
(27 particles/cell/species for the ambient plasma) and 
an electron skin depth, A ce = c/u> ve = 9.6A, where 
Wpe = (47re 2 n e /w e ) 1 / 2 is the electron plasma fre¬ 
quency and A is the grid size [5, 8]. 

The particle number density of the jet is 0.741rib, 
where rib is the density of ambient (background) elec¬ 
trons. The average jet velocities in the two simulations 
are Vj = 0.9798c corresponding to Lorentz factors of 
5 (2.0 MeV). The jets are cold (v? th = 'ej] th = 0.01c) 
in the rest frame of the ambient plasma where c is 
the speed of light. Electron-positron plasmas have 
mass ratio m p /m e = 1. The mass ratio of ion and 
electron is rrii/m e = 20. The electron (positron) 
and ion thermal velocities in the ambient plasmas are 
u® h = = 0.1c and cj h = 0.022c, respectively. The 

time step At = 0.013/tCp e . 

The electrons are deflected by the transverse mag¬ 
netic fields (B x , By) via the Lorentz force: — e(v x B), 
generated by current filaments (J z ), which in turn en¬ 


hance the transverse magnetic fields [5, 6, 8]. The 
complicated filamented structures resulting from the 
Weibel instability have diameters on the order of the 
electron skin depth (A ce = 9.6A). This is in good 
agreement with the prediction of A ~ 2 1 / 4 c7 t 1 1 { 2 /w P e ~ 
1.188A ce = 11.4A (lsj. Here, 7th ~ 1 is a thermal 
Lorentz factor. At the earlier time smaller current 
filaments are generated. However, in the electron- 
positron jets because of larger growth rates, the cur¬ 
rent filaments are coalesced in the transverse direction 
in the nonlinear stage as shown by Nishikawa et al. 
[13, 14]. The longitudinal current (J z ) in the electron- 
positron jet shows significantly more transverse vari¬ 
ation than in the electron-ion jets [13, 14], 

Current filaments resulting from development of the 
Weibel instability behind the jet front are shown in 
Fig. 1 at time t = 28.8/Lv pe for unmagnetized electron- 
positron ((a) and (c)) and electron-ion ((b) and (d)) 
ambient plasmas. Figures la and lb show the y- 
component of the current density in the z — x plane 
at y = 43A. The maximum values of J y are (a) 15.6 
and (b) 23.9, respectively. The arrows show the per¬ 
turbed currents (J z , J x ). Figures lc and Id show the 
^-component of the current density in the x — y plane 
at z = 230A. The arrows show the magnetic fields 
(B x , By) created by J z and their maximum values are 
(c) 54.7 and (d) 95.6, respectively. Based on these 
figures the maximum values of the generated current 
densities (magnetic fields) are larger in the electron- 
ion ambient plasmas. The ambient ions are responsi¬ 
ble for the larger maximum current densities. 

It should be noted that as shown in the previous 
work for the different jet types ((a) electron-positron 
and (b) electron-ion) injected into ambient plasmas 
with the jet composition, the maximum current den¬ 
sities have different structures [13, 14]. Current fila¬ 
ments resulting from development of the Weibel 
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FIG. 3: Velocity distributions of jet electrons two cases as Fig. 1 at t = 28.8 /ui pe ((a)and (b): electron-positron ambient 
plasma, (c) and (d) electron-ion ambient plasma. Jet electrons are binned as a function of , yv± ((a) and (c)) and 7 V 11 
((b) and (d)), and "fv±, where 7 = (1 — (tijf + Vx)/c 2 ) -1 ^ 2 . The blue and red curves show the distributions of injected 
and shocked jet electrons. 


instability behind the jet front (at z = 230A) 
are obtained at time t = 28.8/w pe for two different 
cases with unnragnetized ambient plasmas. The max¬ 
imum values of J z are (a) ±54.7 (as same as Fig. 
la), (b)—123.6, respectively. The electro-positron jets 
show electron (negative) and positron (positive) cur¬ 
rent filamentations, since both species contribute to 
the Weibel instability. On the other hand, at this 
simulation time mainly electron jets generate (nega¬ 
tive) current filamentations. In the electron-ion jet at 
this time only electrons are contributing in exciting 
the Weibel instability, electron current channels are 
dominate, therefore J z has the negative value. 

The differences in the generated current densities 
(J z ) between the different ambient plasmas are seen 
more clearly in the x-component of the generated 
magnetic fields as shown in Fig. 2. The ampli¬ 
tudes of B x in the electron-ion ambient plasma (b) 
are much larger than those in the electron-positron 


ambient plasma (a). The distances between peaks in 
(b) are clearly larger than those in (a). This seems 
to come from the fact that the ambient ions are af¬ 
fecting the growth of the Weibel instability. These 
differences in the amplitudes and structures will af¬ 
fect the emission, which needs to be calculated from 
an ensemble of all particles (electrons and positrons). 
Since magnetic fields are generated by the Weibel in¬ 
stability “jitter” radiation need to be considered by 
the deflected electrons (positrons) [15, 16, 17]. In or¬ 
der to examine the saturated magnetic fields we need 
much longer simulations, which are in progress at the 
present time. 

The Weibel instability is aperiodic, i.e., w rea i ~ 0 
(convective) [12]. Thus, it can be saturated only by 
nonlinear effects and not by kinetic effects, such as 
collisionless damping or resonance broadening. Hence 
the magnetic field can be amplified to very high val¬ 
ues. This characteristic is seen in Fig. 2b. The time 
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evolution of the B x shows that its amplitudes grow as 
the jet propagates without any significant oscillations 
in the frame of the jet. 

Simulation results show that smaller scale current 
filaments appear immediately behind the jet front. 
Smaller filaments merge into larger filaments behind 
the jet front in the non-linear stage. This phenomenon 
is seen in Fig. 2. In particular Figure 2b shows a large 
and long wavelength just behind the jet front. 

It should be noted that a long simulation with 
electron-ion jet with 7 = 15 show that the electron 
Weibel instability switches to the ion Weibel instabil¬ 
ity [9]. The ion current channels with the ion skin 
depth contribute to accelerating electrons [4, 7, 9]. 

The acceleration of electrons has been reported in 
previous work [4-10, 13, 14]. Figure 3 shows that the 
cold jet electrons are accelerated and decelerated. As 
expected, at this time jet electrons in the electron- 
positron ambient plasma (b) are thermalized more 
strongly than those in the electron-ion ambient plasma 
(d). The blue curves in Figs. 3b and 3d are close to the 
initial distribution of injected jet electrons (half of jet 
electrons). We also see that the kinetic energy (paral¬ 
lel velocity U|j ss t>j) of the jet electrons is transferred 
to the perpendicular velocity via the electric and mag¬ 
netic fields generated by the Weibel instability [11, 
12], The strongest transverse acceleration of jet elec¬ 
trons accompanies the strongest deceleration of elec¬ 
tron flow and occurs between z/A = 210 — 240. The 
transverse velocity distribution of jet in the electron- 
positron ambient plasma (at the peak value) is slightly 
more accelerated than that in the electron-ion ambi¬ 
ent plasma as shown in the red curves in Figs. 3a 
and 3c. The strongest acceleration takes place around 
the maximum amplitude of perturbations due to the 
Weibel instability at z/A ~ 230 as seen in Figs la and 
lb. 


III. SUMMARY AND DISCUSSION 

We have performed self-consistent, three- 
dimensional relativistic particle simulations of 
relativistic electron-positron jets propagating into 
unmagnetized electron-positron and electron-ion 
ambient plasmas. The main acceleration of electrons 
takes place in the region behind the shock front [4-5, 
7-10, 13, 14, 17]. Processes in the relativistic colli¬ 
sionless shock are dominated by structures produced 
by the Weibel instability. This instability is excited 
in the downstream region behind the jet head, where 
electron density perturbations lead to the formation 
of current filaments. The nonuniform electric field 
and magnetic field structures associated with these 
current filaments decelerate the jet electrons and 
positrons, while accelerating the ambient electrons 
and positrons, and accelerating (heating) the jet and 
ambient electrons and positrons in the transverse 


direction. 

The effects of the different ambient plasmas have 
been investigated in this report. In spite of the lo¬ 
cal larger current filaments generated in the electron- 
ion ambient plasma, the total jet electron acceleration 
is smaller than that in the electron-positron ambient 
plasma. 

The growth rates depend on the Lorentz factors of 
the jet as suggested by the theory [ 12 ], E-fold time is 
written as r ~ ^/ 7 sh/w pe , where 7 s h is the Lorentz fac¬ 
tor of the shock. The simulation results show that the 
growth time becomes larger with the larger Lorentz 
factor [18]. 

Other simulations with different skin depths and 
plasma frequencies confirm that both simulations have 
enough resolution and the electron Weibel instability 
is characterized by the electron skin depth [ 8 ]. 

An additional simulation in which an electron-ion 
jet is injected into an ambient plasma with perpendic¬ 
ular magnetic field shows magnetic reconnection due 
to the generation of an antiparallel magnetic field gen¬ 
erated by bending of jet electron trajectories, conse¬ 
quently jet electrons are subject to strong non-thermal 
acceleration [ 10 ]. 

The generation of magnetic fields both with and 
without an initial magnetic field suggests that emis¬ 
sion in GRB afterglows and Crab-like pulsar winds 
could be either synchrotron or jitter emission [15, 16, 
17]. The filament sizes appear to be smaller than can 
produce observable variations in intensity structure. 
However, this small size can mean that the deflection 
angle, a ~ eB j_Ab/ 7 m e c 2 , of particles by Weibel fil¬ 
aments is smaller than the radiation beaming angle, 
Ad ~ I /7 [ 12 ]. Here Ab ~ A ce , eB±/m e c < fl e , and 
the ratio S ~ a/A9 < D e /ui pe will be less than one 
when the cyclotron frequency is less than the plasma 
frequency. Thus, when ambient magnetic fields are 
moderate, i.e., the cyclotron frequency is less than the 
plasma frequency and <5 < 1 , the emission may corre¬ 
spond to jitter rather than synchrotron radiation [15, 
16, 17]. 

The fundamental characteristics of relativistic 
shocks are essential for a proper understanding of the 
prompt gamma-ray and afterglow emission in gamma- 
ray bursts, and also to an understanding of the par¬ 
ticle reacceleration processes and emission from the 
shocked regions in relativistic AGN jets. Since the 
shock dynamics is complex and subtle, more compre¬ 
hensive studies are required to better understand the 
acceleration of electrons, the generation of magnetic 
fields and the associated emission. This further study 
will provide insight into basic relativistic collisionless 
shock characteristics needed to provide a firm phys¬ 
ical basis for modeling the emission from shocks in 
relativistic flows. 
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